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Abstract 

We study the flavor physics in two extensions of the quark sector of the Standard 
Model (SM): a four generation model and a model with a single vector-like down- 
type quark (VDQ). In our analysis we take into account the experimental constraints 
from tree-level charged current processes, rare Kaon decay processes, rare B decay 
processes, the Z — > bb decay, K, B and D mass differences, and the CP violating 
parameters |-, ek and a^K- All the constraints are taken at two sigma. We find 
bounds on parameters which can be used to represent the New Physics contributions 
in these models (A^, \ b t f and \ s t f in the four-generation model, and £7^, Uj, s and U s d 
in the VDQ model) due to all the above constraints. In both models the predicted 
ranges for asi (the CP asymmetry in semi-leptonic decays), AMp, B(K + — > 7r + i/p), 
B(Kl — > ifivv) and B(Kl — > [j,£l)sd can be significantly higher than the predictions 
of the SM, while the allowed ranges for a^K and for ArriB s are consistent with the 
SM prediction. 

1 Introduction 

Extensions of the quark sector modify many features of the Standard Model flavor physics. 
In particular, CKM unitarity is violated, and there are new sources of flavor-changing 
neutral-currents (FCNC) and of CP violation. We here analyze the flavor physics of two 
such extensions: a fourth generation of fermions and a single vector-like representation of 
quarks added to the SM. 

These two models share many common aspects. Both models have the same number 
of mixing angles and phases, which can be taken as parameters of a 4 x 4 unitary matrix. 
Both models violate the 3x3 CKM unitarity in a similar manner, and introduce new 
sources of CP violation. Although originating from different Feynmann diagrams, many 
of the expressions for the experimental observables in the two models are practically the 
same (except for numerical differences). The features of CKM non-unitarity and new 
CP violating parameters are of particular interest, since they offer us a behavior which is 
qualitatively different from that of the SM. Also, the two models provide new contributions 
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to FCNC processes. Since FCNC are highly suppressed in the SM, they provide a useful 
tool in searching for New Physics and in constraining it. 

In order to constrain the models, we consider the new contributions to rare K and 
B decays, neutral meson mixing mass differences, CP violating parameters in the K and 
in the B systems, and the process Z — > bb. To obtain the numerical results, we scan 
the allowed ranges of the mixing parameters. For each point in this parameter space, we 
check that all the above mentioned constraints are obeyed. This gives us bounds on the 
possible values of the mixing parameters and the correlations between them. We also use 
the scan to give predictions for various observables in these models, and compare them to 
the predictions in the SM. Finally, we analyze the predictions of these models to various 
processes that are not measured yet, such as K + — > 7i + uu, Kl — > 7r°z/z/, B — > X s £ + £~, 
AM D and AM Bs . 

The paper is organized as follows. In section |2| we describe the analysis for the four 
generation model, and in section [| we describe the analysis for the VDQ model. Each 
of these sections includes some background about the model, descriptions of the various 
constraints used in the analysis, and the details of the numerical results. Finally, in section |4] 
we discuss our conclusions. 

2 A four generation model 
2.1 Background 

Currently there is no known fundamental principle which fixes the number of SM gener- 
ations. Regarding experimental constraints, as of today there is no conclusive evidence 
that excludes a fourth generation. There are, however, two problematic issues. First, the 
invisible decay width of the Z° boson clearly indicates the existence of exactly three light 
neutrinos. The existence of a fourth neutrino with mass m„4, > 45 GeV is not excluded by 
the data, but it requires some mechanism that will give the new neutrino a large mass while 
keeping the masses of the SM neutrinos small. Second, the four-generation model has some 
difficulties in explaining the electroweak precision measurements. Various analyzes of elec- 
troweak precision measurements (e.g. pj and |J) differ in their conclusions regarding the 
implications for a fourth generation. We assume that these measurements are consistent 
with four generations. Specifically, our analysis can be viewed as complementary to the 
one in ref. 0, taking into account also the mixing in the quark sector. 

The constraints we consider are the following: charged-current tree-level decays, the 
branching ratios B(K + — > 7r + b>v), BiK^ — > [ifi)sD and B(B — > X s £ + £~), the mass differ- 
ences ArriK, Am Bd , Am Ba and Am^, the CP violating parameters Sk and a^ K , and the 
partial decay width of Z — > bb. 

An extensive work on this subject can be found in Q. We add to it the B — > ipKs 
constraint, which was not available at the time, and the Z — > bb constraint. Also, we 
treat B{B — > X s £ + £~) in a more careful manner. In addition, we update the experimental 
bounds used in || and we consider the entire range of possible mixing angles and phases. 
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Another work which considered the tree-level decays, a rough analysis of the Z — > bb 
process and the B meson system constraints, can be found in [|J. 

2.2 The model 

We consider a model with an extra generation of chiral fermions added to the SM in the 
simplest possible way. The representations of the fourth generation fermions are identical 
to the representations of the usual three SM generations. In such a model, there are nine 
additional parameters compared to the SM: four new masses, and five new mixing angles 
and phases. The origin of the new mixing angles and phases comes from the fact that the 
CKM matrix is now a 4 x 4 matrix. This means that it consists of nine physical parameters: 
six mixing angles (compared to three in the SM) and three complex phases (compared to 
one in the SM). We use a specific parametrization of the 4x4 CKM matrix 0, in order 
to incorporate all the correlations in the analysis which we perform. The mixing angles 
will be referred to as #12, #13, #23, #14, #24 and #34, and the phases as 0i3, 0u and 024- In the 
limit of vanishing new mixing angles, one gets the usual mixing angles #12, #13, #23 and the 
CP violating phase 0i 3 of the SM. 

We deal only with the quark sector of the model, which includes the SM quarks as well 
as the new quarks (to be referred to as (f b')). The new leptons do not play any role in 
our analysis, and we do not discuss them. In particular, we do not restrict ourselves to 
any specific mechanism that gives rise to the high mass required for the fourth neutrino. 

The effects of the fourth generation enter the processes we consider through loop dia- 
grams: the new quarks t' and b' can now appear in the loops. In most of the constraints, 
the fourth generation contribution enters almost exclusively through the parameters A^, 
Ajf and A^/ d , where 

^ = V kl V km- (1) 

These parameters can be used to evaluate the effects of New Physics contributions. 

2.3 The Constraints 

The analyzes done in concerning electroweak precision measurements with four gen- 
erations indicate that only a restricted range of m t i — my is allowed: 

\m t > ~m b '\< 85 GeV (at 95% CL). (2) 

Another constraint on the masses comes from the combined results of direct measurements 
at CDF and DO collaborations ||, which rule out my < 175 GeV, assuming that the 
FCNC decay mode is dominant. We assume hereafter m t i > my > 175 GeV. Also, since 
a perturbative approach is assumed to be valid, one must require m t > , my < 500 GeV. 
In view of the above constraints, we consider in our analysis two extreme cases of the 
new fermion masses. The first uses the highest possible mass values (in view of the above 
constraints): m t i = 500 GeV, my = 470 GeV. The second uses m t i = 200 GeV, my = 170 
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GeV. For each of the two cases (m t / = 500 GeV or m t i = 200 GeV), any alternative choice 
for my (which obeys the limitation given in eq. (fj)) can cause only minor changes to the 
results we present in this work. The actual calculations for both cases (m t / = 500 GeV and 
mt' = 200 GeV) are basically the same, with only numerical differences. In order not to 
repeat the presentation twice, we give a detailed analysis only for the case m t > = 500 GeV. 
The results for the second case (with m t > = 200 GeV) are briefly summarized in section 



Measured SM tree-level processes are not affected by the existence of a fourth gener- 
ation, and can therefore be used as in the SM to constrain elements of the CKM matrix. 
This enables us to bound the magnitude of the matrix elements V U d, V us , V u b, V c a, V cs and 
V c b. The parametrization of Vckm means that we can translate these constraints into 



Parameter 



Mean value 



Sigma 



Ref. 



\Vud\ 
\Vus\ 

I Kb I 
\V cd \ 

\Vcs\ 

\V*\ 

Rb 

Br(K L — > (ip,)sD 
Br{K+ -> ti+vD) 
Br(K L -> n°vi>) 
B(B -> X s e+e~) 

SK 

S 

Attik 
Am Bd 
Am Bs 

a*pK 

\M? 2 \ 



0.9735 
0.2205 
0.00349 

0.224 

0.996 

0.041 
0.21653 

< 3.75 x 10" 

< 5.07 x 10' 
< 5.9 x 10" 

< 1 x 10~ 5 
2.28 x 10" 3 
17.2 x 10~ 4 
3.489 x 10" 15 GeV 

3.2 x 10~ 13 GeV 
> 9.87 x 10~ 12 GeV 
0.8 

< 6.3 x 10~ 14 GeV 



-t) 

-10 
-7 



0.0008 
0.0018 
0.00076 
0.016 
0.024 
0.002 
0.00069 
95% CL 
95% CL 
95% CL 
90% CL 
0.013 x 10~ 3 
1.8 x 10~ 4 
0.008 x 10- 15 GeV 
0.092 x 10~ 13 GeV 
95% CL 

0.1 
95% CL 



MM 



Table 1: Values of experimental-related input parameters used in the analysis. 



bounds on some of the mixing angles. Taking the constraints as described in Table [I] at 
two sigma, and using unitarity of the 4 x 4 Vckm matrix, we deduce the following ranges: 

sin #12 = 0.22 , 
sin 9 23 = 0.041 ±0.004 , 
sin #i3 = 0.00349 ± 0.00152 , (3) 
< sin #14 < 0.085 , 
< sin #24 < 0.25 . 

These bounds are independent of the new quark masses. The parameter sin# 34 and the 
CP violating phases 0i 3 , i4 and 2 4 remain unconstrained at this stage. 
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We now consider loop processes. These depend on the Inami-Lim functions [19 





X (x) 



Y Q {x) and S (x) defined in 



Experimental inputs are given in table [I] and theoretical 
ones m table g. Most of the expressions for the various constraints in the four generation 
model are obtained in a simple manner from the relevant expressions for the SM, by 
considering additional diagrams with t' quarks in the loop. 



Parameter 


Value 


Ref. 


!s d v 7 B B d 

f Bd y/B Bd 
fK 

B K 


0.23 ± 0.04 GeV 
> 0.2 GeV 

1.15 ±0.06 

0.16 GeV 
0.85 ±0.15 


m 
§ 

1 

[7| 



Table 2: Values of decay constants and bag parameters used in the analysis. 



The rare Kaon decays which we consider are Kl — > lip. and K + — > 7i + uu. A bound on 
B(Kl — > hP)sd (the short-distance contribution to the dispersive part of B(Kl — > fip)) 
can be extracted from the experimental data as described e.g. in ref. 0. The expression 
for this quantity in the four generation model is given by 

B(K L -> fifc)sD = 6.32 x 1(T 3 x [Y NL Re\ s c d + rj t Y Y (x t )Re\f + r] t , Y Y {x t ,)Re\¥] \ (4) 

where Y^l = (3.5 ± 0.6) x 10 -4 |2fJ| represents the charm contribution. In our analysis 
we neglect the QCD correction factors r] t Y and t]t> Y , since they are close to unity. The 
expression for B(K + — > 7r + vv) in the four generation model is given by 

B(K + -> tt+z/z/) = 1.55 x 10- 4 x \\f v f X (x t ) + \ s t f^X{x t ,) + X NL \f f , (5) 

where X^l = (9.8 ± 1.4) x 10~ 4 represents the charm contribution. In this case we again 
neglect the QCD corrections which are close to unity. We also neglect the uncertainty in 
Xnl- 

Next we discuss mass differences in the various neutral meson systems. The four gen- 
eration expression for the short-distance contribution to Am^- is given by 

Am sn = GHfem^ ^^ + x f v K So{xt) + X f v K So (x t ,) + 

+ 2\f\f v %S (x c ,x t ) + 2\?X?r&S (x„x 1 ,) + 2\fA s /v%S (x u x t ,)]. (6) 

Assuming that the long-distance contributions can be at most of the order of the experi- 
mental value Am^ p , we take this constraint as < Am)p <7x 10~ 15 GeV. 

The contributions to the short-distance part of the D° — D° meson mixing amplitude 
M\2 in the SM are very small compared to the current experimental bound (by several 
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orders of magnitude). The leading diagram for Mf<f in the four generation case contains 
two b' quarks in the loop, and is given by 

I M SD I - G F m wfp B D m D it r * y* |2„ c („ \ ( 7 \ 
M 12 I - 77) 2 \ V cb' V ub'\ Vb'Oo{X b ') ■ (7 J 

Since the long-distance contributions to M 12 are also estimated to be small compared to 
the experimental bound, we demand in the four-generation case < 6.3 x 10~ 14 GeV, 

where the experimental bound was taken as in table [I] with rjy = 0.56. The experimental 
bound we use is weaker than the one usually considered (e.g. in since it does not assume 
5 = and uses instead the available experimental data ||18|| . From the parametrization of 
Vckm, \V u b'V c b'\ 2 ~ sin 2 ^14 sin 2 ^24, so eq. (0) gives the bound 



sin #14 sin 6> 2 4 < 



4.1 x KT 3 m t , = 500 GeV 
9.1 x 10- 3 m t > = 200 GeV 



The four generation expression for Am# d is given by 
G 2 F rn 2 w f 2 B B Bd m Bi 



Am 



B d 



6tt 2 



xfr ] fS Q {x t ) + X h f^S {x t> )+2Xf\fr ] f tl S G {x u x v ) . (9) 



We take the QCD factors as rjf = 0.55, r]f t , = 0.5 and r$ = 0.54. 

The expression for ArriB s is the same as that for Am,B d , the only change is to replace 
d — > s in all places. Using the input of tables ffl and 0, one obtains 



Xf V ?S (x t ) + Xf^SoM + 2X\ s X\?r 1 ZS Q {x u x tl ) 



> 19.6 . (10) 



X b t d2 ^S (x t ) + X b /r,*So(x t ,) + 2X bd X bd r ) ^ l S (^ x t 
Next we consider the Ek constraint. In a four generation model, one has 

£k = \c F B K x Im[Xf *\f S Q (x c ) + Xf\fS (x t ) + X s f\$S Q {x t ,) + 

+ 2X sd *Xf \l S (x c , x t ) + 2Xfxfr,«S (xe, a*) + 2Xfxfr,«S (x t , a*)] . (11) 

where C e = Gj,.Tr^f^m K /{&r 2 V2AmK) = 3.8 x 10 4 . For the SM QCD corrections we 
take J7J r)J* = 0.57, = 1.38 and rj%. = 0.47. For the new QCD corrections we take |3| 
r§ v = 0.57, t$ = 0.6 and rfe = 0.5. 

Another constraint, a^x, comes from the decays B — ► ^fC. In the SM, this quantity 
measures to an excellent approximation sin 2(3, where (3 = aig(— X bd / X bd ) is one of the 
angles in the unitarity triangle. In the presence of New Physics, this is modified ]22[] to 
sin(2/3 — 28 d), where (3 is defined as in the SM and 8d is the New Physics phase of the 
mixing amplitude M 12 . In the case of four generations, 8& is given by 

28- ,rJ Ml2 \- anrfl I V*S (x t ,)X b f + 2 V *S (x t ,x t ,)Xyxf 
26d " [mW) ~ " arg { 1 + ^ (x t )A- 2 ' ' (12) 
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where M 12 and M^ M are the mixing amplitudes for the B meson system in the four- 
generation model and the SM respectively. For the experimental data we use the world- 
average of a^K-, given in table [l], at two sigma. Recently, a new preliminary result regarding 
the B — > ipK measurement was published by the BABAR collaboration p3 |. As this is 
still only a preliminary result, we did not add it to our analysis. However, these results are 
not expected to have a dramatic impact on the results obtained here. 

In the following subsections we discuss in more detail two additional constraints used 
in our analysis, those coming from T(Z — > bb) and from B{B — > X s £ + £~). 



2.4 TiZ 



bb) 



In this section we follow the derivation and the notations of ref. [23], and add to it the 
contributions of the fourth generation. The decay rate for Z —>■ qq can be written as 



T(Z - qq) 



a 



Ifis 2 r 2 



-rn z (\a q \ 2 +\v q \ 2 )(l + 5 g V)(l + 5« QED )x 



x (1 + 5% CD ){1 + 5«)(1 + 8l QCD ){\ + 6 b ) . (13) 



Here a is the 



\ and Cyy are, respectively, the sine and cosine 



We use a caret for quantities given in the MS renormalization scheme, 
electromagnetic fine structure constant, s 
squared of the Weinberg angle, a q = 2/| and v q = 2/| — AQqS^ are the relevant axial and 
vector coupling constants respectively, and the 5 terms are corrections due to various high 
order loops. In the calculation of Rb, most of the 5 terms cancel out and are therefore 
irrelevant for our discussion. The terms which remain are S^q CD , which consists of QCD 
contributions to the axial part of the decay, 5^ which consists of kinematical effects of 
external fermion masses (including mass-dependent QCD corrections), and 5b which is 
different from zero only for q = b and is due to the Z —>■ bb vertex loop corrections. 

The contributions to the corrections 5>1q CD originate from doublets with large mass 
splitting, and are given in the four-generation case by 



5 q 

u tQCD 



a.. 



— [atf(fh) + <h'f{lH') + avfbn,)] 




f(fi f ) « log | - J - 3.083 + 



0.346 0.211 



4 ' 



(14) 



where /i 2 = 4m 2 /m%. The corrections to the Z — > bb vertex in the four generation case 
are given approximately (for m t i = 500 GeV) by 



5 h 



10 



-2 



— 0.0154|T4 fo | 2 - 0.148|V t / 6 | 2 



5| + - 2 )'" ; " 12 



(15) 



The corrections 5 q are given in ref. |24], and are not affected by the fourth generation. 



Using the expressions given in ref. |24j for R s , R C} R U} together with eqs. Ql4]) and (0), one 
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obtains 



R b = (1 + 2/R s + 1/R C + 1/R U ) 



1 + 



3.5837 1 1 



1 — 0.01541 Kbl 2 — 0.148|l^/ fc | 2 



(16) 



Eq. ( |16"D together with the experimental bound taken from table p] at 95% CL imply 

2 i 1 2 

|Vt&| + 9. 5 1 Vt'b\ < 1.14. This can be used, together with the unitarity of the third column 
of Vckm and the parametrization of Vckm, to get 



. . . f 0.14 m t / = 500 GeV 

8111034 ~( 0.3 m t , = 200 GeV (17) 

This bound becomes weaker as m t > is lowered, and is only very weakly dependent on 
my. In the degenerate case, m t = m t >, one cannot put any limitation on the mixing angle 
#34, as expected. 

This constraint contributes significantly to the analysis. It excludes the maximal mixing 
solution found in ||, as well as some of the other mixing solutions suggested there. In 
addition, it probably implies that the up quark mass matrix has some special structure. 
This can be seen if we recall the assumption m t < m t / < 500 GeV, which means that 
there is no hierarchy between the third and the fourth generation up quark masses. One 
then naively expects large mixing between the third and the fourth generation, so that 
sin #34 = 0(1). However, the bound on sin #34 (eq. ([[7])) implies that in the four-generation 
model this naive expectation is not fulfilled, and the mass matrix has some non trivial 
structure. 



2.5 B(B -> X s i + i~) 

In previous analyzes (e.g. ||), only the leading contributions of Z-mediated diagrams to this 
branching ratio were considered. However in light of the improvement in the experimental 
data, we use a more detailed analysis of this quantity. In this section we use mainly the 
derivation and the notations of [pOfl . In our analysis, however, we completely neglect QCD 
corrections. The decay can be described at low energies by using the following effective 
Lagrangian : 



C eff = —j=X? C 2 (sc) v _ A (cb) v _ A + C 7 ^m b sa^(l + 75)^+ 



V2 



•57T 



+C 9V (sb)y_ A {ll) v + C 10A {sb)y_ A {ll) A ] 



where (qiq2)v±A = 917^(1 ^ 75)^2 and m& is the b quark mass. The coefficients C 7l) C^v 
and Ciqa originate in the SM from electroweak magnetic penguin diagrams, Z and 7 pen- 
guin diagrams and W box diagrams. The coefficient C2 contributes to the decay via an 
intermediate cc loop. In order to get the total branching ratio, one should integrate the 
differential branching ratio given in [20| over the entire allowed range, which we take as 
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4m? 



< s < 1 (here s = qf+^/mf, where q is the total dilepton momentum). The resulting 



u 

bound (for electrons in the final state), when neglecting QCD corrections, is 



6s 12 



w, 



cb\ 



{(0.5|(7 9 | 2 + 0.5|C 10 | 2 + 0.93|C 2 | 2 + l.lRe[C*C 9 ] + 125.6|C 77 | 2 + 



+4Re 



4.3Re[C 2 C* 7j ] < 42.4, (19) 



where Cg = —Cgy and Cio 



2tt 
a 



Cioa- 



The values of the coefficients C 2 , C 77 , and C 10 are given in the SM by 

C 2 = l, 



Cg = - In 
9 



C 



10 



m B 
Y (x t ) 
sin 2 W ' 



(20) 



sin 2 fl 



IV 



where D' (x t ), Y (x t ), Z (x t ) are Inami-Lim functions. When considering a model with a 
fourth generation, the coefficients C 2 , C 77 , C 9 and Ci change as follows: 



C 2 = 1 + 



A t 6 *' 



C 77 (M W , / 



1 / A bs \ 

~2 + j& D ' oM ) 
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In 



M 



m B 



1 ^/ 

1 + aF. 



^bs 

Y (x t ) + -£Y (x t 



a* 



sin 2 6w 



C 



\bs 

Y (x t ) + -T^Y (x t 



10 



sin 2 # 



YV 



(21) 



We use in our analysis the bound of eq. (|19"D, with the coefficients as in eq. (^T 



2.6 Numerical results 

In the numerical analysis we take all the constraints that were described in the previous 
subsections, and scan over the allowed region in the parameter space (namely, the mixing 
angles according to bounds given in eq. ([|) and ([H]) and the three phases between and 
2tt). The scan is performed by randomly choosing the values of the mixing parameters 
(using a uniform distribution), and checking for each such point in the parameter space 
that all the constraints are met at two sigma. When extracting from the scan the allowed 
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ranges for various quantities (such as mixing angles and predictions for observables) we 
used only 95% of the data points given by the scan results. In this section we analyze the 
results for the case of m t i = 500 GeV. The results for m t i = 200 GeV are given in the 
following section. 

According to the scan, some of the new mixing angles and phases have restricted ranges 
(see figure |l|): 

7T 

< <Pi3 < 2 ' 
0<sin# 14 < 0.031 , 

0<sin# 24 < 0.033 , (22) 
0<sin# 3 4<0.14. 

The remaining new phases (cfiu and ^24) can take the entire scanned range. 




Figure 1: Histograms of (a) 13 , (b) sin# 14 and (c) sin# 2 4 for 60,000 data points, for 
m t , = 500 GeV and m v = 470 GeV. 



Note that in the SM the allowed range (at two sigma) for ^13 is tt/5 < ^13 < An/9. 
This is quite close to the range of this phase in the four generation model. But while in the 
SM the allowed range is dominated by the Ek constraint together with Am^ and Amg s , 
in the four-generation model it is influenced by a correlated effect of several constraints 
(and not only the three that dominate in the SM case). 

In order to see in a clearer way the impact of the various constraints on New Physics, 
one can look at the allowed regions of the parameters A*, d , \\f and X t f. Their scatter plots 
are given in Figure |2|. Our analysis shows that there is very little correlation between 
the various constraints. This means that each of the parameters is effectively constrained 
only by a small number of constraints. We now consider each of the parameters and the 
constraints that influence it most. 

1. \ s t f: The bound on ImX^ comes mainly from the B(K + — > 7i + uu) constraint. The 
real part, ReXff, is mainly constrained by B(K L — > ^lJi) sd . The general shape of Xp 
is mainly due to the Ek constraint. 
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Re^ d xio 4 Re^ d ReX* 



(a) (b) (c) 

Figure 2: Scatter plot of the New Physics parameters (a) A|, d , (b) \ b t t and (c) \ b t f with all 
the constraints, for m t i = 500 GeV and my = 470 GeV. 



2. Aj, : |Aj/ 1 is mainly constrained by the Am^ constraint. The a^x constraint is 
responsible for the area excluded in the upper-right and lower-left regions in figure 
|2(b)| . We also examined how an improvement in the experimental error for a^x will 
affect the influence of this constraint. As it turns out, the excluded area is almost 
insensitive to the size of the error, and even an error as low as 0.01 (10 percent of 
the current error) will have very little effect in the (Re\\f, ImX^f) plane. A change 
in the central value from 0.8 to 0.7 will have some influence on the excluded area. In 
general, the higher the central value is, the larger is the excluded area. However, the 
modifications in the excluded area are in any case not drastic. This means that the 
coming experimental results for a^x are not expected to change the bounds on the 
fourth-generation flavor parameters in a significant manner. 

3. A^f: This parameter is constrained mainly by the B — > X s £ + £~ bound, together with 
the tree-level constraints and the Z — > bb constraint. It is little influenced by the 
ArriB a bound. Our bounds for this parameter are better than those previously quoted 
in the literature, since we use an updated (and stronger) experimental constraint for 
B -> X s £+£-. 

A good way to compare the four-generation model to the SM is to look at the predictions 
given by the two models for various quantities. The results are summarized in table [| We 
give here several notes regarding these predictions: 

1. The predictions for |M^|, B(K~l — > fj,p,)sD, B(K + — > 7t + vu) and B(K~l — > t[°uu) in 
the four-generation model can be significantly higher then the SM predictions. In 
the cases of B(Kl — > fJ>p)sD an d B(K + — > tt + vv) they can span the entire range 
implied from the experimental constraints. 

2. The CP asymmetry in semi-leptonic decays is approximately given by the model- 
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B(K L — > A*A*)sd 


B(K L -> tt u w) 




B(B -> X s l+l~) 


experiment 


< 3.75 x 10" 9 


< 5.9 x 10~ 7 


< 5.07 x 10~ 10 


< 1 x 10~ 5 






[i 7 ^ fil v 1 n~ n 


[4 7 11 nl v 1 n -11 
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[15,32] 




[-0.47, +0.28] 





Table 3: Predictions for different quantities in the SM and in the four-generation model 
compared to the experimental bounds, taken at two sigma. The category of 'favored range' in 
the four-generation model refers to the range that contains 95% of the data points. 



independent expression [30] 



sin 2^ (23) 



(r 12 /M 12 ) 5M \M l2 /M™\ 

The experimental bound on this quantity is obtained by calculating the world-average 
of clsl = (0.2 ± 1.4) x 10~ 2 from [[2T)|-|2"9"|] , and taking into account the theoretical 



predictions for (T 12 /M 12 ) 1 « -(0.79±0.27) x 10~ 2 (by an updated scan according to 
the data in ref. P0"|). Although our scan improves the bounds which were obtained in 
[§], the four-generation model prediction for this quantity can still be about an order 
of magnitude higher than that of the SM. Also, this quantity in the four generation 
model may have a different sign from that predicted by the SM. The four-generation 
prediction for clsl is quite far from the current experimental bound. In case the 
experimental bound improves by about a factor of 3 for the upper bound and 4.5 for 
the lower bound, this bound will become significant in the analysis. 

3. Taking all constraints at two sigma, the a^x prediction in the four-generation model 
covers the entire range allowed by the B — > ipK constraint. The SM prediction covers 
the range 0.6 < a^K — 0.95. However, the difference between the SM and the four- 
generation model regarding this quantity can be clearly seen if we perform the scan 
without including the B — > ipK constraint. Then the four-generation model allows 
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the range —1 < a^ K < 1, while the SM allowed range is only 0.4 < a^ K < 0.95. Still, 
from the latest experimental results it is clear that the value of cl^k lies at the higher 
part of the allowed range, so this information is of limited impact. 

4. The Ai7iB s prediction in the four-generation model is similar to the prediction of 
the SM. This means that detection of Am^ outside the SM range will be a major 
problem not only for the SM but also for the four generation model. 

2.7 The case of m t > = 200 GeV, m v = 170 GeV 

In this case, the bounds that are obtained for the mixing parameters are significantly weaker 
than for the case of m t i = 500 GeV, my = 470 GeV. The fact that the strongest effects 
are obtained for the heavier mass of t' can be explained as follows. Roughly speaking, the 
diagrams that we consider put bounds on terms of the form X^y^x-f), where yi(x t ') is 
some Inami-Lim function. Since yi(x t i) grows with x t > = m t / 2 /mw 2 , the heavier the mass 
is, the stronger is the bound on the parameter X 1 ™. 

The histograms for the phase 0i3 and the angles sin # 14 and sin #2 4 are qualitatively the 
same as for the case m t > = 500. The resulting ranges of the new mixing angles and phases 
are given by: 

< 013 < ^ ' 

O<sin0 14 < 0.051 , 

0<sin# 24 < 0.056, (24) 
0<sin# 3 4<0.3. 

The remaining new phases (0i 4 and ^24) can again take the entire scanned range. 

The bounds on the New Physics parameters X^f, X^f and X s t f are presented in figure 
|3]. In this case, X\f is constrained mainly by the tree-level decay processes, Z — > bb and 
unitarity, and not by B — > X s £ + £~ as in the case of m t > = 500 GeV. Other than that, 
one can see that the bounds for the cases of m t t = 500 GeV and m t i = 200 GeV are 
qualitatively the same, with obvious numerical differences. 

The predictions for various quantities in the case of % = 200 GeV, my = 170 GeV are 
quite close to the those in the case of m t ' = 500 GeV, my = 470 (table Basically this 
means that the allowed ranges of the four-generation model for these quantities is almost 
independent of the new quark masses. 
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Figure 3: Scatter plot of the New Physics parameters (a) \\?, (b) A^f and (c) Xff with all 
the constraints, for m t i = 200 GeV and my = 170 GeV. 



3 Vector— like Down— type Quarks (VDQ) 
3.1 Background 

We consider a model in which a single VDQ is added to the SM. Vector-like quarks 
transform as triplets under the SU(3) C symmetry, and both their left-handed and right- 
handed components transform as singlets under the SU(2) L symmetry. VDQs are predicted 
by various extensions of the SM, such as grand unified theories based on the Eq Lie algebra. 
Also, as explained in section III models with VDQ exhibit interesting features such as the 
violation of CKM unitarity and the related appearance of FCNC contributions at tree-level, 
and the appearance of additional CP violating phases. 

The constraints we consider in our analysis are the following: charged-current tree- 
level decays, the branching ratios B(B — > X s £ + £~), B(K + — > 7T + uu) and B(Ki — > fifi)sD, 
the mass differences Am B[j and Am Bs , the CP violating parameters Ek, — and a^x, and 
precision electroweak measurements related to the Zbb coupling. The scan procedure that 
we use is similar to that used for the four generation model. 



This model was previously studied in the literature (see [31] and references therein). 
Regarding Rb and A b FB (the forward-backward asymmetry of the b quark) in this model, 
see [Q. Recent works similar to the one presented here were performed by [33]. We 
update the experimental bounds used there, taking into account the new results from Belle 
and BABAR collaborations regarding the B — > ipK measurements. We also present the 
predictions of this model for various observables and compare them to the predictions of 
the SM. While this paper was in final stages of writing, another article 0] was published 
on this subject. 
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3.2 The model 



We consider a single VDQ, denoted b', added to the SM. This means that the mass matrix 
in the down sector is now 4x4, and it is diagonalized by a 4 x 4 matrix V. The down sector 
charged-current interactions depend on the 3x4 upper submatrix of V, which plays the 
role of the CKM matrix, but are otherwise unchanged. The neutral-current interactions 
are given in this model by 

£nc = -^--Z^ li Ul 1 - Wl,U a pD^ - 2 sin 2 9 W J* M ) , (25) 

Z COS o\y 

where a sum over repeated indices is implied (a, (3 — d, s,b), and we used the definitions 




In eq. ([25|), is the electromagnetic current, which contains also terms with the new 
quark b' . Note that eq. (|25"D contains FCNC In all the processes we consider, the leading 
New Physics contributions come from the tree-level FCNC that appear in eq. ( |2li| ) through 
the quantities U a p- These New Physics contributions usually compete with contributions 
coming from SM loop processes. Other New Physics effects, due to b' quarks in loop 
diagrams, are naturally highly suppressed compared to the tree-level contributions, and 
can be safely neglected. We can thus view this model at low energies as having the same 
particle content and interactions as the SM (completely ignoring the extra b' quark), but 
having a non unitary CKM matrix. 

Note that the quantities Ub s , Um and U s d play two important roles. First, they indicate 
the amount by which the three-generation CKM matrix deviates from unitarity. Second, 
they represent the strength of the new contributions to FCNC. These quantities play the 
same role as — X^, in the four-generation model. 

The matrix V is not a general unitary matrix; some of the phases in it can be removed 
by change of basis. It can be parameterized by nine parameters: six mixing angles and 
three phases. All these parameters appear also in the upper 3x4 submatrix. As in the 
four-generation model, we use the specific parameterization of || for the 4x4 matrix V, 
in order to incorporate all the correlations in our analysis. The mixing angles are again 
referred to as 9u, #13, #23, #14, #24 and #34, and the phases as ^13, (fiu and 024- 



3.3 The constraints 

The constraints obtained from the tree-level decays are the same as in the four-generation 
model, described in section |2]3|. The expression for B(B — > X s £ + £~) in the VDQ model is 
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given by eq. (|I9[) , with the coefficients C 2 , Cg and C w taken as 



r - 1 ^ 
62 " 1_ aF' 



* 8 (M w \( U bs \ Yo(x t ) ( ( 1 \ C/ 6s / 



c 



Y (x t ) + C U2Zx 



10 



sin 2 6 



The coefficient C*7 7 remains unchanged. 

The expressions for the remaining constraints are very similar to those used in the four 
generation case. They are taken as in P5|| , but with the input parameters given in tables 
HI and 0. For the — constraint we use only the ranges of the parameters given in [35|. 



3.4 Numerical results 

The numerical scan is performed in the same way as in the four generation case. The 
results of the scan show that the CP violating phases and 024 can be in the entire 
scanned range (from to 2tt). Yet, as in the four generation case, the phase 0i3 has a 
restricted range, < 0x3 ^ 71 ■ This results agrees with [33|. A histogram for this phase 
can be seen in figure |4(a)| . The mixing angles #14 and 824 also have restricted ranges: 

0<sin# 14 < 0.011 , 

0<sin# 24 < 0.011 . (28) 
The restricted regions are due to combinations of several constraints. For example, the 




Figure 4: Histograms of (a) 13 , (b) sin #14 and (c) sin 824 for 75,000 data points in the 
VDQ model 



restricted range of sin 8 24i is mainly due to the correlated effect of the K L — > /j,/! and j 
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constraints. The remaining phases, 024 and 0i4 can be in all the scanned range. The 
mixing angle #34 can be in the range < sin #34 < 0.12 (see [|33|j). 

We now examine the allowed regions of the parameters U sd , Ubd and Ub s , which represent 
the New Physics contributions in this model. Their scatter plots are given in Figure |5]. 
They are in very good agreement with the results of |33|, except for changes which are 
related to the recent a^x and B — > X s £ + £~ measurements. As in the case of the four 
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Figure 5: Scatter plot of the New Physics parameters (a) U s d, (b) Um and (c) Ub s with all 
the constraints. 



generations, also in this case each of the parameters is effectively constrained only by a 
small number of constraints. We now consider each of the parameters: 

1. U s d- The bound on ImU sd comes mainly from the j constraint. The real part, ReU s d, 
is mainly constrained by B(Kl — > ^Jj)sd- 

2. Ubd'- \Ubd\ is mainly constrained by the Am^ constraint. The a^x constraint is 
responsible for the area excluded in the upper-right and lower-left regions in figure 
|5(b)| . We examined how an improvement in the experimental error for a^x will affect 
this constraint. As in the four generation case, the results show that the excluded 
area is almost insensitive to the size of the error, and there is little sensitivity also 
to the central value. This means that the coming experimental results for a^x are 
not expected to change the bounds on the flavor parameters of the VDQ model in a 
significant manner. 

3. Ubs- This parameter is mainly constrained by the B — ► X s £ + £~ bound. It is little 
influenced even by Am^. The bounds we obtain are better than those obtained in 
previous works due to the improved experimental bound (see table [[]) that we 



use. 



The last step in comparing this model to the SM is to consider the predictions given by 
the two models for various quantities (see table §). We give here several notes regarding 
these predictions: 
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1. The predictions for B(K L — > fip,)sD, B(K + — > -n + vu) and B(K L — > 7r°z/z/) in the 
VDQ model can be substantially higher than the predictions of the SM. The predic- 
tion for B{K L — * lip) sd can actually span the entire range allowed by the experi- 
mental constraints. 

2. The prediction for B(B — > X s £ + £~) in the VDQ model can reach the current exper- 
imental bound. This is in contrast to the prediction in the case of low a^x values, 
as presented in | 33| . In a similar way, the predictions for BiB — > Xd£ + £~) are also 



large in this model. 

The scan results for the CP asymmetry in semi-leptonic decays (as given in table 
|3|) improve the bounds which were obtained in Note that the prediction in the 
VDQ model may have a different sign than the SM prediction. In order for the 
asL constraint to have a significant impact on the parameters of the VDQ model, 
the experimental results must improve by about an order of magnitude. Such an 
improvement is not expected in the near future. 

The behavior of the predictions for cl^k are similar to the four generation case. When 
taking all constraints at two sigma, the results for the SM and for the VDQ model are 
practically the same. But when performing the scan without including the B — > i/jK 
constraint, the VDQ model allows the entire range for a^K, while the SM allows only 
a very restricted range. Still, from the latest experimental results it is clear that the 
value of a^K lies at the higher part of the allowed range, so this information is of 
limited interest. 

As in the four generation case, the Am^ prediction in the VDQ model is similar 
to the prediction of the SM. This again means that detection of Amg 8 outside the 
predictions of the SM will be problematic also in the VDQ model. 



4 Discussion and Conclusions 

We considered constraints from tree-level decays, electroweak precision measurements, the 
decay Z — > bb, rare K and B decays, CP violating parameters in the K and in the B 
systems and mass differences in various neutral meson systems in order to obtain bounds 
and predictions for the four generation model and for the VDQ model. The constraint on 
the four generation model from Z — >■ 66 proves to be very powerful, as it excludes a signif- 
icant portion of the parameter space. Improvement in the B(B — > X s £ + £~) measurement 
provides better bounds on the New Physics parameters in both models. The new experi- 
mental data for a^x also affects our results, causing the analysis for the VDQ model to be 
quite different from the results of fl33f . Furthermore, according to our analysis additional 



improvement of this measurement is expected to have only minor effects on both the four 
generation and the VDQ models. Improvement of the bounds for these models can thus be 
expected mainly by better determination of other measurements such as B(B — > X s £ + £~), 
B(K + — > 7i + uu), asL, and Amg s . Of course, both of the models can never be strictly 
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excluded by flavor physics, since by proper adjustment of the new parameters one can 
always reduce them to the SM. In the four-generation model, this can be done by taking 
the mixing between the fourth generation and the three SM generations to be very small 
(sin 6*34, sin Ou, sin 6* 24 — > 0). The result is decoupling of the fourth generation from the 
other three. In the VDQ model, this can be achieved by taking the mass of the new quark 
to be very high. Since the new mixing angles result from the diagonalization of the mass 
matrix, they are roughly given by 6 i4 m so by taking ^> rrii we can make all the 
new mixing angles vanish. 

We find that the two models which we consider share many common features. They 
both contain new sources of CP violation that arise from a 4 x 4 CKM-like matrix, and 
both predict for various observables values that can be higher than those predicted by the 
SM. Another aspect that the two models share is that they both introduce new FCNC 
contributions. In both models these new FCNC contributions are naturally small, though 
they can still be larger than the SM contributions. This was shown to lead to substantial 
increase in the values the two models predict for various quantities, compared to the SM 
predictions. Despite these (and other) similarities between the four-generation and the 
VDQ models, there are also significant differences between them. One of the differences 
is related to the fact that the mechanism that introduces the new FCNC sources is not 
the same in the two models. The new FCNC contributions in the four-generation model 
come from loop diagrams with the new fermions in the loop. These diagrams are roughly 

proportional to (^j^j ~ O(10). In the VDQ model, on the other hand, the leading new 
FCNC contributions come from tree-level diagrams, with a coefficient of Cjj2Z ~ 0(100) 
(see e.g. ||33||). Thus the new contributions for the FCNC processes in the two models differ 
by an order of magnitude. 

This order-of-magnitude difference leads to various effects. One of these regards the 
unitarity relation X bd + X bd + X\ d + X\t = in the four-generation model, or X bd + X bd + X\ d - 
Ubd = in the VDQ model. As was already discussed in the literature (see e.g. [|] and 
references therein), in both cases one gets unitarity quadrangles instead of the unitarity 
triangle that exist in the SM. When examining the results of the numerical scan (see table 
f|), it is clear that the possible shapes that the unitarity quadrangle can take is very different 
in the two models. The parameter r is defined by 



One can see that in the four-generation model, X d can be significantly larger than the 
corresponding SM quantity X bd (up to a factor of 3.5 for m t i = 500 GeV and a factor of 70 
for m t i = 200 GeV), while in the VDQ model Um is only allowed to be about 15% of X bd . 
This difference can be traced back to the order-of-magnitude difference between the new 
FCNC contributions in the two models. As a result of this, also the ranges of the angle 
j3 are significantly different in the two models. While in the VDQ model the unitarity 
quadrangle can be only slightly modified compared to the SM triangle (due to the limited 




four generation model 



VDQ model. 



(29) 
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Model 


(3 (degrees) 


r 




SM 




22 - 33 





four- 


generations (me 


= 500GeV) 


f) - 57 281-360 


< 3.5 


four- 


generations (me 


= 200GeV) 


- 92, 230-360 


^ 70 




VDQ 




2 - 38 


^ 0.16 



Table 4: Allowed ranges of various parameters related to the unitarity quadrangle, when all 
constraints are taken at two sigma. r is defined in eq. (|29|). 



size of r and in the four-generation model the shape of the unitarity quadrangle can 
be completely different than that of the SM unitarity triangle (especially in the m t i = 200 
GeV case). 
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